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Abstract 

Traditional few-body descriptions of pion production use integral equations 
to sum the strong interactions nonperturbatively. Although much physics 
is thereby included, there has not been a practical way of incorporating the 
constraints of chiral symmetry into such approaches. Thus the traditional few- 
body descriptions fail to reflect the underlying theory of strong interactions, 
QCD, which is largely chirally symmetric. In addition, the lack of chiral 
symmetry in the few-body approaches means that their predictions of pion 
production are in principle not consistent with the partial conservation of 
axial current (PCAC), a fact that has especially large consequences at low 
energies. We discuss how the recent introduction of the "gauging of equations 
method" can be used to include PCAC into traditional few-body descriptions 
and thereby solve this long standing problem 

I. INTRODUCTION 

Approximate chiral symmetry is an important property of quantum chromodynamics 
(QCD) and should therefore be an attribute of any effective description of strong inter- 
action processes. Yet most nonperturbative effective descriptions of pion production in 
few-body processes (NN — > nNN, ttN — > ttttN, 'yN — > ttN, etc.) have been based on 
traditional few-body approaches which are not consistent with chiral symmetry - even if 
they are based on chiral Lagrangians (i.e. the input functions are chirally invariant). The 
essential idea of a traditional few-body approach is to take as input two-body t matrices 
and pion production vertices, and then use integral equations to sum the multiple-scattering 
series nonperturbatively. A feature of such an approach is that the input can be constructed 
phenomenologically without the need to specify an explicit Lagrangian for the strong inter- 
actions. Despite the large amount of physics taken into account (e.g. all possible pair-like 
interactions are included), the lack of approximate chiral symmetry in this approach results 
in a pion production amplitude that does not obey PCAC. Although this is a particularly 
serious problem at low energies where low-energy theorems apply, the absence of PCAC is 
undesirable at any energy because of its inconsistency with QCD. In this contribution we 
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show how the gauging of equations method |],|2| can be used to construct a nonperturba- 
tive few-body approach that is consistent with chiral symmetry and whose pion production 
amplitude obeys PCAC. 



II. TRADITIONAL FEW-BODY MODEL OF NN -» ttNN WITHOUT PCAC 

We shall base our discussion on the example of the relativistic nNN system for which a 
traditional few-body approach has recently been developed [§,|J. For simplicity of presen- 
tation we shall treat the nucleons as distinguishable. The few-body approach to the nNN 
system provides a nonperturbative simultaneous description of the processes NN — > nNN, 
NN — > NN, and nNN — > ttNN within the context of relativistic quantum field theory. In this 
case the integral equations are four-dimensional and can be expressed symbolically as^ 

T = V + Vg t T (1) 

where, for distinguishable nucleons, T, V, and Q t are 4x4 matrices written as 

t — I Tnn T n \ ( V NN T \ r _ ( D \ , v 

T -[t n Ty V -\T Gfiy Gt -{0 G w°G )- (2) 

The elements of matrix T are defined as follows. T is a 3 x 3 matrix whose elements are 
Alt-Grassberger-Sandhas (AGS) amplitudes (generalised to four dimensions) describing the 
process nNN — > nNN. Note that the following "subsystem-spectator" labelling convention 
is used: A = 1 or 2 labels the channel where nucleon A forms a subsystem with the pion, the 
other nucleon being a spectator, while A = 3 labels the channel where the two nucleons form 
the subsystem with the pion being the spectator. In a similar way, Tjviv is the amplitude for 
NN scattering while T N and T N are 3x1 and 1x3 matrices whose elements T AA r and T Nfl 
describe AW — > nNN and nNN — > NN, respectively. For simplicity of presentation we shall 
neglect connected diagrams that are simultaneously NN- and nNN- irreducible. Then the 
elements making up the kernel matrix V specified in Eq. (0) take the following form: 

V m = V^ E - A (3) 

where V^ E is the nucleon-nucleon one pion exchange potential and A is a subtraction term 
that eliminates overcounting. T is a 3 x 1 matrix with 

2 

^A = E^-fi (4) 

1=1 

where Fi = fidj 1 consists of /j, the vertex for A, — * nNi, and dj, the Feynman propagator 
of nucleon j ^ i. The subtraction term B in Eq. (ffl) likewise eliminates overcounting. T 



x In this paper it should be understood that all amplitudes and currents with external nucleon 
legs are actually operators in Dirac space and need to be sandwiched between appropriate Dirac 
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FIG. 1. Diagrams making up the subtraction terms, (a) W W7r , (b) W^n, (c) Wnn, (d) X, (e) Y~, 
and (f) £>. The dark circles represent the following two-body amplitudes: (a) full irir t-matrix, (b) 
one-nucleon irreducible ttN t-matrix, and (c) full NN t-matrix minus the NN one-pion-exchange 
potential. 



is the 1x3 matrix that is the time reversed version of J 7 , Gq is the ttNN propagator, and 
I is the 3x3 matrix whose (A,/z)'th element is 8\ tlt . Finally the propagator term Q t is 
a diagonal matrix consisting of the NN propagator D , and the 3x3 diagonal matrix w° 
whose diagonal elements are tid^ , £2^iA an d t^d^ 1 , with t\ being the two-body t matrix 
for the subsystem particles in channel A (for A = 1 or 2, t\ is defined to be the 7rN t matrix 
with the nucleon pole term removed). The subtraction terms A and B are defined with the 
help of Fig. |l] as follows: 

A = W n7T + W^ N + W NN + X + Y'- B'G B' (5) 

where W^ N = W wN + PW^ N P, Y' = Y + PYP, and B' = B + PBP, P being the nucleon 
exchange operator. 

Despite the rich amount of physics incorporated into the NN — > ttNN amplitude within 
this model, it becomes evident from the following discussion that this amplitude cannot 
satisfy PCAC. 



III. NEW FEW-BODY MODEL OF NN -» ttNN WITH PCAC 

A. Choosing appropriate degrees of freedom 

Although we do not specify the exact Lagrangian behind our new few-body approach, 
we do assume that this underlying Lagrangian is chirally invariant (up to a small explicit 
chiral symmetry breaking term). In turn, the chiral invariance of the Lagrangian puts a 
strong constraint on the nature of the fields that can be used to construct a practical few- 
body model. For example, if one would like to follow the ttNN model above and have only 
pions and nucleons as the degrees of freedom, then the underlying chiral Lagrangian would 
necessarily involve an isovector pion field £ that transforms into functions of itself under a 
chiral transformation. Unfortunately, this can happen only in a nonlinear way |J: 

(6) 

where 9 is the vector of three (infinitesimally small) rotation angles. Thus one pion trans- 
forms into two pions, two pions transform into four pions, etc., a situation that would make 
it difficult to formulate a chirally invariant few-body description (i.e. a description whose 
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exposed states have a restricted number of pions - in our case or 1). Alternatively, we 
can follow the example of the linear sigma model and choose an underlying Lagrangian that 
involves, in addition to pions and nucleons, the isoscalar field a. In this case the a, together 
with the pion field 0, transform under the chiral transformation in a linear way: 

$^$+6a, o^o-e-4>. (7) 

Thus the four-component field <fi = (a, <p) transforms into itself - a situation that is ideal 
for a few-body approach. We shall therefore adopt the latter approach and treat pions and 
sigma particles on an equal footing. For the few-body nNN model of Sec. 2, this means 
mostly a formal change where the usual (isospin) three-component pion is replaced by a 
four- component one. There is, however, one new aspect in that terms with a three-meson 
vertex (7nra) now need to be included. For example, one will need the new subtraction term 
illustrated in Fig. |. Nevertheless, the equations of Sec. 2 retain their structure, the only 
essential change being an increase in the size of the matrices like T, T/v and T/v to take into 
account the introduction of a aNN channel. With these modifications we obtain few-body 
ttNN equations that are consistent with an underlying chiral Lagrangian. 
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FIG. 2. New subtraction term in the few-body irNN model where pions and sigmas are treated 
on an equal footing. 



B. Coupling the axial vector field everywhere 

Now that we have achieved consistency with an underlying strong interaction chiral 
Lagrangian, our next goal is to construct a few-body nNN model whose axial current is 
partially conserved. The problem at hand is analogous to the one of constructing a conserved 
electromagnetic current for a few-body system whose strong interactions are consistent with 
an underlying Lagrangian that conserves charge. In this case it is well known that exact 
current conservation is obtained by coupling an external electromagnetic field to all possible 
places in the strong interaction model. In a similar way, one would expect to obtain a 
partially conserved axial current by coupling an external axial vector field to all possible 
places in the strong interaction model. Until recently, however, what has not been known is 
the way to achieve this complete coupling for a few-body system whose strong interactions 
are described nonperturbatively by integral equations. Fortunately, the solution to this 
problem presented recently in the context of electromagnetic interactions [0,0, is based on 
a topological argument and therefore applies equally well to the present case of an axial 
vector field. 

The basic idea is to add a vector index (indicating an external axial vector field) to all 
possible terms in the integral equations describing the strong interaction model. Thus, in the 
case of the ttNN system where the structure of the integral equations is given as in Eq. ([[]), 
coupling an external axial isovector field gives 
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(8) 



which can easily be solved to give a closed expression for T M : 

r» = (i + Tg t )v»(i + g t T) + Tg?T. (9) 

T M is a matrix of transition amplitudes Tfc N , Tjy A , T^ d , etc. (for the moment, we suppress 
the isovector index in these amplitudes). Here we are particularly interested in the transition 
amplitude Tfc N as it is closely related to the pion production amplitude we seek. 

It is easy to see that the transition amplitude Tfc N given by Eq. (pj) has the axial vector 
field attached everywhere in Tjvjv except on the external nucleon legs M. Including these 
external leg contributions then gives the complete axial vector transition current for NN — > 
NN: 



f NN = (r^x + r%d 2 )T NN + iwfar? + d 2 r%) + t» n (10) 

where Tf is the axial vertex function of nucleon i. The input to the irNN equations con- 
sists of the two-body t matrices t i7 the pion production vertices /j, and the single particle 
propagators d\. Thus the input to Eq. (|) also includes the axial transition currents if, /f , 
and axial vertex function Tf . In order for these input quantities to be consistent with an 
external axial vector field being attached everywhere, they must be constructed to satisfy 
the Axial Ward-Takahashi (AWT) identities |J. This is easily achieved by restricting the 
form of the corresponding bare quantities |L|. Thus, for example, needs to satisfy the 
AWT identity 

<Z M ir(M = i [d-\k) l5 + Kdr\p)] t a - iUmlf?(k,pK(q) (n) 

where p and k are the initial and final momenta of the nucleon, q = k — p, f^(k,p) is the 
irNN vertex for a pion of isospin component a, t a is an isospin 1/2 matrix, is the mass 
of the pion, and f w is the pion decay constant. With the other input quantities constructed 
to satisfy similar AWT identities, it can be shown that the two-nucleon axial current given 
by Eq. fllCf ) satisfies the AWT corresponding to exact PCAC: 

QlijNN(klk 2 ,PlP2) 

= i [(l5t a )iTNN(h - q, k 2 ;pip 2 ) + T NN (k 1 k 2 ;p 1 + q , p 2 ) (j 5 t a ) i] 
+ i [{ r Y5t a ) 2 T NN (k 1 , k 2 - q;pxp 2 ) + T NN (k 1 k 2 ;p 1 ,p 2 + q){ r Y5t a ) 2 ] 
- lUmlT^ihk^p^d^q) (12) 

where T^ is the amplitude for ttNN — > AW and in its time reversed form, is just the the 
pion production amplitude that we are seeking. Note that the axial current contains a 
pion pole : 

j%N( k ik 2 ,Pip 2 ) = j%N( k ik 2 , Pip 2 ) + 4(g)i^(?)^ (^i^2,PiP2), (13) 

where j^ N has no pion pole and F£ is the pion decay vertex function. Thus an alternative 
way of obtaining the pion production amplitude Tjv , is to take the residue of Eq. (|T0D at 
the pion pole. 
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C. Physical content of the new model 



Not only does the pion production amplitude T/vo obey exact PCAC, but it contains a 
very rich amount of physics that goes beyond what is included in the traditional few-body 
approach of Sec. 2. A few of the infinite number of new contributions are illustrated in 



FIG. 3. Examples of new pion production mechanisms included in the few-body model with 
PCAC. 

Finally we would like to stress the practicality of our approach: PCAC and inclusion of 
an infinite number of new pion production mechanisms has been achieved through a simple 
closed expression, Eq. (||), involving scattering amplitudes T obtained from a traditional 
few-body model, Eq. ([I]). Our method also does not depend on the model used for kernel 
V - rather than Eq. (Q), we could equally have used the simpler case of an AW one-pion- 
exchange potential. 
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